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In this work we use functional integral techniques to examine the nearest neighbour attractive 
Hubbard model on a quasi-2D lattice. It is a simple phenomenological model for the high-Tc 
^ ' cuprates that allows both extended (non-local) s- and d-wave singlet superconductivity as well as 

\ mixed symmetry states. The Hartree-Gor'kov mean field theory of the model has a finite tempera- 

ture phase diagram which shows a transition from pure s-wave to pure d-wave superconductivity, via 
a mixed symmetry s + id state, as a function of doping. Including Gaussian fluctuations we exam- 
ine the crossover from weak-coupling BCS superconductivity to the strong-coupling Bose-Einstein 
$Ih ' condensation of composite s- or d-wave bosons and comment on the origin and symmetry of the 

O ' pseudogap. 
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I. INTRODUCTION 



I [ The physics of the high-Tc cuprates is unusual and fascinating for a number of reasons. In this work we explore 
'"O two of them: the symmetry of the superconducting pairing state and the existence of the pseudogap state between 
the insulator and the superconductor. 
^ , In recent years the experimental evidence has become overwhelming that the cuprate high-Tc superconductors 
possess an anisotropic d-wave pairing state Angle resolved photoemission (ARPES) experiments ||,^ and mea- 
surements of the temperature dependence of the penetration depth |^-^ both strongly support this scenario, but 
perhaps the most convincing evidence comes from the elegant measurements of Wollman et al. |^] and of Tsuei et 
al. since they do not depend on the microscopic physics of the energy gap, but instead depend only on the order 
parameter phase. 

' Whilst there is a consensus on the d-wave nature of the dominant part of the order parameter, there is continuing 
^\ , controversy over the possible existence of a sub-dominant s-wave contribution P,p^. The existence, or otherwise, of a 
' subdominant s-wave component has profound consequences for the nature of the microscopic pairing mechanism. For 
On , example, antiferromagnetic spin fluctuations lead to attraction in the d^2_y2 pairing channel, but are pair breaking 
in the s-wave channel . Similarly the Hubbard model with a positive (repulsive) on-site interaction U may have a 



dx2_y2 paired ground state, but would presumably not support s-wave Cooper pairs. In the case of YBa2Cu307 the 
orthorhombic crystal symmetry makes some non-zero s-wave component inevitable, but a large s-wave component 
would be difficult to reconcile with either of these pairing mechanisms. On the other hand, pairing mechanisms based 
. on electron-phonon interactions, polarons, or other non-magnetic excitations (e.g. excitons, acoustic plasmons) could 
(— I ■ be compatible with either s-wave or cZ-wave pairing states [ p^ . In these models, whether s-wave, d-wave or mixed pairs 
Q , are more strongly favoured would depend on details of the model parameters, and could even vary from compound to 



O ■ compound. Indeed there is some evidence that the n-type cuprate superconductors are s-wave |13 (or at least they 
^ ' have no zeros in the gap |A(k)| on the Fermi surface). This would imply that either the pairing mechanism is different 

■ '~j ' for the n- and p-type materials, or that the mechanism allows both s-wave or d-wave ground states depending on the 

f\ band filling. 

^ . The second remarkable feature of the cuprates is the "pseudogap" occuring in underdoped materials. Transport 
■ - - ' experiments @, ARPES |,|], NMR and optical conductivity all detect a loss of spectral weight at low 

energies in the one-particle excitation spectrum below a certain temperature, T* > Tc. Furthermore, the k-specific 
nature of the ARPES technique has shown that the pseudogap has nodes at the same places as the superconducting 
gap. These data and the fact that the coherence length is known to be very small in the cuprates (^(0) ~ 10 — 20 A 
in the Cu02 planes) points to a strong-coupling scenario in which superconducting fluctuations play a major role. 

Although there is as yet no concensus on the exact nature of the fluctuations, it is generally believed that the 
opening of the pseudogap at T* is due to the formation of an incoherent pairing state above Tc. 

Two competing scenarios have emerged to explain this phenomenon. In the first, exemplified by the work of Emery 
and Kivelson [h9|, the relevant fluctuations are in the phase. The stiffness to the phase fluctuations is ns(0)/m* 
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(where 71^(0) is the zero temperature superfluid density and m* is the effective quasiparticle mass), and where this is 
smah, such as in the underdoped cuprates, the true Tc, where phase coherence is estabhshed, can be weh below the 
mean field transition temperature which is identified with T*. 

The second scenario is best illustrated by the works of Nozicres and Schmitt-Rink ||2^ (which predates the discovery 
of the high-Tc compounds), Dreschler and Zwerger and Sa de Melo, Randeria and Engelbrecht (see Ref. p^ ] 
for a review) . By considering both magnitude and phase fluctuations in s-wave superconductors in the vicinity of the 
critical temperature, these authors examined the crossover from BCS-like superconductivity of weakly bound pairs 
to the Bose condensation of tightly bound pairs as U, the strength of the interaction, was varied. The results are 
summarised in Fig. ^ For weak pairing the fluctuating and mean field theories are essentially the same. However 
as U is increased large differences become apparent: whilst the mean field critical temperature, T^^, keeps getting 
larger, the effects of fluctuations cause the true T^ to flatten off, for the continuum model, or even to decrease in the 
case of the lattice model. The difference between T*^^ and Tc gives an estimate for the regime where fluctuations 
dominate and may correspond to the pseudogap regime in which bosonic pairs exist but are not condensed into a 
superfluid. 

In this paper we will address both the pairing symmetry and the fluctuation issues in a systematic way. We 
broadly follow the functional integral approach taken by Dreschler and Zwerger |^ and Sa de Melo, Randeria and 
Engelbrecht , but explicitly including c?-wave pairing. In section |l| we introduce the attractive nearest neighbour 
Hubbard model as a simple phenomenological model of the cuprates and derive fro m it a Landau-Ginzburg- Wilson 



(LGW) functional describing s- and d-wave superconducting fields. In section HI we examine the BSC-like mean 



field theory of the model by finding the saddle point solutions of the LGW functional which show s- and c?-wave and 



mixed symmetry phases. Finally, in section IV, we include Gaussian fluctuations about the saddle point in order to 
study the strong coupling bosonic limit. This limit is especially interesting in this model because of the competition 
between preformed pairs of different symmetries. We also comment on the symmetry and origin of the pseudogap. 

II. LANDAU-GINZBURG- WILSON THEORY 

To investigate the symmetry of the order parameter and the nature of the bosonic limit for anisotropic pairs, we 
consider a simple phenomenological model of the cuprates, the attractive nearest neighbour Hubbard model on a 2D 
square lattice. It is the simplest model that allows s-wave, d-wave or mixed symmetry pairing states to occur. Its 
Hamiltonian is, 

H = Ho + Hi, (1) 

where Hq contains the chemical potential and a nearest neighbour hopping term. The interaction Hamiltonian 
describes an electron-electron attraction on nearest neighbour sites: 

HI = -VJ2n^n„ (2) 

{ij) 

where denotes summation over all of the bonds between the nearest neighbour sites. For the sake of simplicity, 

we do not include an on-site interaction term. 

In order to investigate the thermodynamics of this model it is sufficient to consider the partition function. In the 
coherent state functional integral representation it is, 

Z = J V[c\c]e^, (3) 



r/3 

S = I dr 







where c and are the Grassmannian counterparts to the Fermi operators c and . 

Unfortunately, the interaction part of the Hamiltonian is quartic in the Grassmann (fermion) fields so the functional 
integral cannot be evaluated exactly. However, by using the Hubbard- Stratonovic transformation ||2^,|2^ the quartic 
problem of interacting electrons can be converted to the more tractable quadratic problem of non-interacting electrons 
coupled to a random (Gaussian weighted) Bose field. In order to make this transformation, the electronic Hamiltonian 
must first be decomposed into a product of fcrmionic bilinears. This can be done in a number of different ways, but 
the purpose of this paper is to study superconductivity so we choose what might be called the Gor'kov decomposition 
in terms of pairing operators: 
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{ij) 



y^tr(4^,,), (4) 



with the pairing matrix written as, 



where cr is the vector of Pauh matrices, and 



i^, = :^(i3,,-^.B,,), (5) 



= (cjTCii - cjj.c^t) , (6) 
7f (cjic,| +Cjtc,|), (7) 

The scalars, are even under both time reversal and parity whilst the vectors, B^j, are odd under both. In the 
manner of Volovik and Gor'kov [2^ , wc interpret Bij and as the annihilation operators for singlet and triplet 
Cooper pairs on the bonds (ij). In terms of these operators the interaction Hamiltonian becomes: 



Ej = -V ^[BlB,,+^\^-^,,), (8) 

<u> 

where we have used tr((Ti) ~ and XxioiGj) — 2Sij. 

Although it would be interesting to study triplet superconductivity, especially in light of the exciting recent discovery 
prf of p-wave pairing in strontium ruthenate, Sr2Ru04, the focus of this work is on anisotropic singlet pairing as seen 
in the high-Tc cuprates. Therefore, for the sake of simplicity, we restrict our attention to the singlet fields. 

Having decomposed the interaction Hamiltonian into a product of fermionic bilinears we make the Hubbard- 
Stratonovic transformation. Tensor valued Bose fields, A^j, are introduced conjugate to the fermion bilinears, Bij. 
(Had we kept the triplet operators we would have needed to introduce an additional vector of triplet pairing fields. 
Ay .) In terms of these the partition function becomes, 

Z = J I?[A*,A] e-^, (9) 
5 = -V^''dT^|A,,(T)|2+ln J 2?[c^c]e^^ 



with the linearised fermions contained in, 

^4(r) {-dr)c,{T)-Ho{c\c) 



rl3 

Sf = / dr 





~F^(A*^.(r)i3,,(r)+A,,(r)4(r) 

(y> 

The Hubbard-Stratonovic fields, A^, are defined on each bond (ij) of the lattice. However, it is generally more 
convenient to work with site-centred combinations possessing a definite symmetry. On a square lattice, every site, i, 
has two fields associated with it, one in the x-direction and one in the y-direction. (We could associate four fields 
with each point but that introduces problems with double counting.) Rather than talk about x- and j/-bond fields, it 
is natural to form linear combinations of them: 



A, = -(A, + A,), (10) 
Ad = ^(A,-A,), (11) 



where Ao(ri) = Au+a for a = x,y according to the direction of the bond {ij). We call these the s- and d-wave 
fields because, in the spatially homogeneous limit, they become the order parameters for (extended) s- and d-wave 
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superconductivity. Unlike the x and y fields, the s and d order parameters transform as irreducible representations 
of the tetragonal crystallographic point group, D^h- In the notation of Ref. they belong to the Aig and Big 
representations of D4^h- The d-wave representation, -Big, corresponds to an order parameter with — symmetry 
as seen in the cuprate materials. 

Proceeding with the calculation, the action is diagonalised by the introduction of Fourier transformed fields: 



c{k) 
Aa(g) 



(12) 
(13) 



in which fc = (ici;„,k) [q — (ia;,y,q)] where the fermion [boson] Matsubara frequencies are w„ — (2n+l)^ — (2j^)^]. 

It is then a straightforward matter to integrate out the vestigial fcrmionic degrees of freedom. The purely bosonic 
action then reads 



5 — 5o + Sh, 

where the non-interacting and interacting terms are respectively, 

5o = -Trln(Go), 

Sb = -^Y. + Trln(l - VGo). 



(14) 

(15) 
(16) 



The first term, 5o, corresponding to the free electron action depends only on the Nambu Green's function matrix. 



Go(fc,A:') = 



:^_( Go(fc) 









^G5(fc) 



(17) 



in which Go{k) — {iujn — ^k) ^ is the free fermion propagator, and ■fk = Ck ^ For nearest neighbour hopping the 
band-energy is ek = -'^tj^i=x.y cos(fci). 

The second term, Si,, which describes the effects of electron-electron interaction, depends not only on the electron 
Green's function but also on the bosonic off-diagonal potential matrix, 



V(fc,fc') = 



V 







A„(fc-fc') 



where the Einstein summation convention has been used for repeated Greek indices and. 



In the limit that k = k' these become, 



C"(k) = Ck.k = 



cos(fca;) -t- cos(fcj,) a 
cos(fca;) — cos(fcj,) a 



a = s.d. 



s, 
d. 



(18) 



(19) 



(20) 



which reflect the point group symmetries of the order parameters. 

Expanding the logarithm in Sh as a power series to fourth order, only the second and fourth order terms are 
non-zero. They are shown diagrammatically in Fig. ^ 

After a trivial rescaling of A by ^ the bosonic action to second order reads. 



(2: 



(21) 



where the particle-particle susceptibilities are, 

xf{q) = +^E^o(^)^o(fc + 'z)Ck.k+qC, 





k+q,k' 



(22) 



= --E 



n/(Ck+a) +n/(^k-a) - 1 



;k+i 
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/■a /-/S 

<^k-H,k+aCk4 



1 k-S J 
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where nf{x) = [e^^ + 1] -'^ is the Fermi function. 

Regardless of the flavour indices, Xo^il) is a maximum at zero frequency and zero momentum. Keeping only 
relevant terms, expansion of the Gaussian action about q = Q gives. 



S ' = — 

2 ^ ' 2m^^ 



The frequency dependence has been kept, to lowest order in uji,, in order to look at the effects of quantimr fluctuations. 
The constant term is given by, 



rO'P = (^SO'P _ ^ ^ tanh 03^/2) 



(24) 



where the weighted densities of states (WDOS) are 

^"''(^) = ^ E r (k)C^(k)<5(e - e^). (25) 

k 

For the nearest neighbour hopping square lattice the WDOS are plotted in Fig. ^. As one can see from the figure the 
d-wave WDOS is strongly peaked around the normal state van Hove peak at half-filling, whilst the (extended) s-wave 
WDOS is suppressed at that point but large near the band edges. 

Note that the square symmetry of the lattice forbids the existence of cross-terms such as r"'^ so that N'^'^ is zero. 
Whilst ^^(k) is even under the exchange of x- and y- momenta, ^'(k) is odd under this operation. The products ('^('^ 
and (^''■(^'^ are both even but C'*C'^ is odd. Provided et is even under x ^ y, iV*'' must be zero. However, if the a; <-> y 
symmetry in is broken, say by an orthorhombic distortion, then r""^ need not be zero. 

The masses come from the second derivatives (with respect to momentum) of the static susceptibilty, Xo'^il) = 
Xo^ {iijJu = 0, q) at q = 0. Generally, for any given band structure, ek, they are. 



2m'ip 8iV 



^^ltel)_!VM^V"(k)C''(k) (26) 



,(k) V a ^? 



k 



+ «,(k)«,(k)(^!^)r(k)C''(k) 
where a prime indicates differentiation and the band masses and quasiparticlc velocities are, 



"'^'^^ vd^)-(?^)- (27) 



m*^(k) \dkidkj J ' \dk 
For the nearest neighbour square lattice only the diagonal terms, m*j(k) and (k), are non-zero and even: 

— ^ = 2tcos{ki), v^{kf = At^sin'^ih). (28) 

Unlike r"'^, the mass in Eq. (E3|) is not purely diagonal in the flavour indices. In addition to the diagonal terms there 
are gradient mixing terms ]29| , |30| , 

^•<"(^)^-"'^^*'(^)^-'"' 

where m'^^ — —fn^g^. They are allowed because they have the symmetry of the square. Under an exchange x ^ y 
both the gradients and the d-wave order parameter change sign, leaving the whole expression unchanged. 

For weak coupling, such terms have little effect except in the vicinity of inhomogeneities such as impurities or 
interfaces where the magnitude of the order parameter is expected to be diminished. In these regions, a gradient in 
the dominant order parameter would drive the subdominant order parameter to be non-zero. So, for instance, near 
the surface of a bulk d-wave superconductor an induced s-wave component might appear. When we consider the 
strong coupling limit we will see the importance of these gradient mixing terms even in the bulk. 
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We now have expressions for r"^ and m^^, leaving only d"^ to be determined. The derivation of (P^ requires 
some care . The Matsubara frequencies are discrete, therefore to expand in them we must first make the analytic 
continuation x°'^{i^vi q = 0) ^ X°'^{'^ + *0+, q = 0) and then expand in u). We find, 

d-^ = \Y. '-^^^^P^N'^^e) - '^N-^{,). (30) 

e 

Like r"^ , the tensor d"^ is diagonal in the flavour indices. The real part of d"^ represents the free propagation of pairs, 
whilst the imaginary part represents the decay of pairs into the continuum of fermionic states. The significance of the 
two terms ^ and their domains of applicability will be discussed when we consider the effects of Gaussian fluctuations 
in section pv| . 

The quartic part of the action is shown in Fig. ^(ii). We associate a susceptibility with the fcrmion loop: 

xf'^' = -74EiG"(*^"'k)|4r(k)c^(k)c(k)c^(k), (31) 



k 

2^E - ^^^^) r(k)c^(k)c(k)c^(k). 



where a, (3, 7, S denote the flavour indices at each vertex. Only the diagrams with an even number of s and d indices 
contribute. Thus, including the rescaling of A by , the bosonic action to fourth order reads, 

5(4) =-iTr(VGo)\ (32) 

= -Us'^A*{qi)A^{q2)Al{q3)A^{qi) x Siqi - q2 + qa - Qa) 
{?} 

-MdE ^d(«i)^d('?2)Ad(93)Arf((74) X S{qi - 92 + 93 - 94) 

{9} 

-2u^'^A*{qi)A^{q2)A*a{q3)Aa{qi) x S{qi - q2 + qs - Qi) 
{9} 

--W E ^s(<?i)Ad(g2)A:(g3)Arf(g4) x 6{qi - q2 + qa - 94) 

E AS(gi)A,(g2)AS(g3)A,(q4) x Siq, - 92 + 93 - 94). 
{9} 



The three coefficients are. 



«. = -^xf- (33) 
". = -^X4--, (34) 
«. = -^xr-, (35) 

where we note that Xi"^'^ = xl'''"'- 

This completes our derivation of the Landau-Ginzburg- Wilson action. In the following sections we will examine the 
saddle point (mean field) and Gaussian fluctuation approximations to this action in order to derive the phase diagram 
of the model and to study the crossover to strong-coupling superconductivity. 



III. MEAN FIELD (LANDAU) THEORY 

In the limit that the superconducting fields are static and spatially uniform, the action, Sb, is related to the Landau 
free energy for interacting s- and d-wave superconductivity: 
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dT = -S^^ = ^T,\^,\^ + i^ri\^d\^ + us\^st + ua\^at (36) 

+2u,|A,nAd|2 + ^ (A,A,A^A^ + c.c.) . 

This form for the free energy is identical to that found by Feder and KaUin |2|] and Joynt who derived it from 
symmetry principles. 

Rationalising the two cross-terms into one by writing the order parameters in terms of their phases and moduli, so 
that 

2u,|A,f|Arf|2 + ^- (A,A,A* A* + c.c.) = 2u^ (^1 + ^ cos(2A0)j lAsHA^p, (37) 

we see that the free energy, JF, depends on the phase difference, A0, between A^ and A^^. 

The saddle points of the free energy, JF, which for a uniform system correspond to Hartree-Gor'kov mean field theory, 
are found by simultaneously minimising T with respect to the magnitudes of the fields and the phase difference, i.e., 

^^^0. (38) 



9|A,| a|Ad| dM 

As a first approximation to solving these equations, we uncouple the two order parameters and, for given V ^ solve 
for Tc as a function of /i for each field independently of the other. In this approximation, Eq. ( p8[ ) gives Stoner-like 
criteria for each of the critical temperatures: 

= 0, a = s, d. (39) 

It is not possible to find analytic solutions to these equations except in the large V limit where we expect Tc to be 
large too. In this case we find, 

lim - 1 - V TV"" (e) , (40) 

e 

which implies that, in this limit, the critical temperatures for both s- and d-wave superconductivity are the same: 

Tc = \- (41) 

Computationally, however, the solution of Tq, = is a trivial matter for all V . The results are shown in Figs. ||(a)-(c) 
for three different values. The maximum value of the ordinate has been chosen as F/4 for each one so that the 
approach to the strong-coupling limit can be followed. 

For small coupling, as in Fig. ||(a), the s-wave solutions have their highest Tc at the edges of the band whilst 
the d-wave solutions have theirs at the band centre. These places correspond to where the respective WDOS of 
each is the largest. The effect of the van Hove singularity in the centre of the d-wave WDOS is to produce a high 
transition temperature in that channel, just as would be predicted by conventional BCS theory in which it is the 
combination yiV""(0) that determines Tc. In fact, it has been credibly suggested by some authors [^l|-p3[ that it is 
the van Hove singularity, observed in most first-principles calculations |^,^, combined with a conventional weak- 
coupling electron-phonon or spin-fluctuation mechanism, that is responsible for the high transition temperatures in the 
cuprate materials. For the extended s-wave channel there is no van Hove singularity and the transition temperatures 
are typically a factor of five lower than the d-wave ones. For small coupling, neither the s- nor the d-wave critical 
temperatures are near the limiting value of V j 4 

As V is increased. Fig. ||(b) , the s-wave solutions begin to occupy the centre of the band and the d-wave solutions 
spread outwards towards the band edges. The maximum Tc's increasingly tend to the strong-coupling value, V jA,. As 
Y increases further towards the strong-coupling limit, Fig. ||(c), the s- and d-wave solutions start to converge. This 
is because the exact shape of the WDOS becomes less important as more of it falls within fc^Tc of the Fermi surface. 
Note that for large V both the s- and d-wave solutions have spread beyond the edges of the band. This, and the 
crossover to local pairs, will be discussed when we consider the effects of Gaussian fiuctuations in the next section. 

Now that the phase diagram of the uncoupled order parameters is known, we turn our attention to the coupled 
case. From Figs. B(b) and (c) we see that there are regions where both the uncoupled s- and d-wave order parameters 
are simultaneously below their critical temperatures. What is the nature of the pairing in these regions? According 
to the values of the model parameters, two cases may be distinguished: (i) For UsUd < (l + ^ cos(2A6')) the 
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pure s- and d-wave phases are separated by a single first order transition which meets the two second order metal- 
superconductor lines at a bicritical point; (ii) For UgUd > (l + ^ cos(2A6')) an intermediate phase of mixed 
symmetry (A^ ^ 0, Ad ^ 0) is permissible, separated from each of the two pure phases by a line of continuous second 
order transitions. These meet the metal-superconductor lines at a tetracritical point. 

In the mixed phase the value of the phase difference, A9, is of paramount importance. It is easily found from Eq. 
(^). Depending on the sign of Ux it can assume either of two values: 

rvr/2 ^.>0 ^ s±^d, 

\Oux<0^s±d. ^ ' 

The two solutions represent fundamentally different phases. For instance, the s ± d phase is invariant under time 
reversal, whilst the s ±id phase breaks this symmetry (time reversal involves complex conjugation). Several authors 
have suggested that the existence of a time reversal symmetry broken order parameter might be associated with spon- 
taneous generation of currents at twin boundaries [p6| or surfaces |3^ . This effect has been observed experimentally 
as a splitting of the zero bias conductance peak in zero magnetic field [Q. No such effect would be observed in the 
s ± d phase. 

Assuming that we are below the critical temperatures of both the s- and d-wave order parameters, all possible 
states of the model can be parametrised by just two dimensionless variables, 



Ud 



_ ux{l + ^ cos(2A6i)) 

Ud 



(43) 
(44) 



where we remember that AO depends on the sign of Ux- 

As the values of these two parameters vary there are transitions between the pure s- and d-wave states and the two 
different mixed phases (separated by a line of first order transitions at Ux — 0)- The stable superconducting phases 
are shown in Fig. ^ as a function of Us and Ux- 

For V = 5t, the variables (Us, Ux) were evaluated along a path in (/x,T) space shown in Fig. ^(b) as a dotted line. 
The resulting values are plotted as the dashed trajectory in ^ Starting in the pure d-wave sector at T = T^, we remain 
in the pure d state even though we are below T^: the s-wave is suppressed by the dominant d-wave order parameter. 
Eventually we cross the second order phase boundary (d — > s ± id) into the time reversal symmetry breaking mixed 
s + id state. With increasing fi the proportion of s-wave in the admixture rises and the d-wave portion is diminished 
until the, now dominant, s-wave totally supresses the d-wave and the second transition [s ± id —^ s) occurs. The 
trajectory terminates in the s-wave sector when T=TJ?. We observe qualitatively identical behaviour for all V. 

The other, s -I- d, mixed phase does not occur spontaneously for any parameters derived from the Hubbard model 



because of the form of the quartic susceptibility, Eq. (31), which ensures that Ux is always positive. To see an 



s -|- d phase it is necessary to include an orthorhombic distortion in the lattice, so that the hoppings in the x and y 
directions are different. Orthorhombic anisotropy, as present in YBCO, introduces an extra term into the free energy 
proportional to A* Ad + A^Ag. With a term like this a finite value for one of the components drives the other to be 
non-zero too. This effect has been studied elsewhere [|36|l so we do not consider it here. 



IV. GAUSSIAN FLUCTUATION THEORY 

It is clear from Figs. H(a)-(c) that, as the coupling V increases, Tc starts to become finite even below the bottom 
(top) of the electronic band at /i = — 4t (-l-4i). For such values of the chemical potential, fermions can only exist 
as bound pairs, not as single particles. A similar phenomenon has been noted for the local s-wave case in Jellium 
by several authors pO|- p^ , p^ who concluded that, in the low density limit, Bose-Einstein condensation of these pairs 
constitutes the mechanism for superconductivity. 

As V increases, the mean field approximation becomes increasingly poor and the effects of fluctuations must be 
included. Specifically we want to know how Tc, as a function of filling, is affected. As it stands, the functional, 
Sb, is too complicated. We therefore employ the Gaussian approximation in which terms of OdAj^) and above are 
neglected. The free energy can be evaluated exactly in this approximation. The action is. 
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There are two important limits to consider: weak coupling/high density (|^| < At) and strong coupling/low density 
> At). Weak coupling corresponds to BCS-like Cooper-pair superconductivity. In this limit the imaginary part of 
d is dominant so the dynamics are relaxational: pairs have a finite lifetime after which they decay into the continuum 
of fermionic states. As the coupling strength is increased, or the density is lowered, the chemical potential can fall 
below the bottom of the band. In this case the imaginary part of d is zero, see Fig. ^, so the pairs become infinitely 
long lived and the physics is that of propagating pairs. This is the Bosc limit. The manner in which the crossover 
comes about is quite different for the s- and d-wave superconductors because of form of the WDOS is different for 
each. The s-wave WDOS (see Fig. ||) shows a sharp band-edge so that the transition from just inside to just outside 
the band is sharp too, whilst the d-wave WDOS (see Fig. ^ goes smoothly to zero so the transition is smooth. 

Although our theory provides an interpolation scheme between the strongly and weakly interacting regimes, its 
validity in the crossover region, near = 4, is doubtful. In this regime, fluctuation corrections pO| of order kpap 
(where kp is the Fermi vector and Oi? is the fermion scattering length), which are vanishingly small in the extreme 
limits, become large signaling a breakdown of the theory. In the preceeding sections we studied the weak limit in 
some detail. In this section we will go on to examine the strong coupling limit, omitting the intermediate regime. 
This is regrettable, especially as it is just this regime that is most relevant to the cuprates, but unavoidable within 
the current formulation. 



As a test of the performance of the Gaussian theory we will first consider the on-site, local s-wave superconducting 
Hubbard model. The Gaussian action for it is, 



where is the local s-wave field and the parameters, r, m and d, are found from those of the nearest neighbour 
extended Hubbard model by replacing V with U and setting (^"^'^^ '*(k) ~ 1. Scaling the fields, we can set the frequency 
coefficient to unity, in which case the action resembles that of free bosons: 



with /ifc = —r I d the bosonic chemical potential, the zero of which signals the onset of Bose-Einstcin condensation, and 
M = md the boson mass. It is well known that as the on-site coupling, U, increases, so too does the mass of the 
bosons. This can be understood by considering the mechanism of the bosonic hopping. The hopping of a pair consists 
of three stages: (i) A pair occupies a single site, (ii) Virtual ionisation - one electron hops onto a neighbouring site, 
gaining a kinetic energy, t, but sacrificing the pair-binding energy, U. (iii) Pair recombination - the second electron 
hops onto the neighbouring site too, gaining a kinetic energy, t, and pair-binding energy, U. The matrix element for 
the whole hopping process, proportional to one over the boson mass, is t^ jU . We now show that our Gaussian theory 
can produce this effect. 

In the extreme large U limit, where the density is very low and U /Tc 3> 1, the electronic chemical potential equals 
half the binding energy of a Bose pair: 



A. On-site (local) pairing 




(46) 




(47) 
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(48) 



In the same limit we find the frequency coefficient. 




(49) 



and the unsealed mass 




(50) 




(51) 
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where the first term in the expansion of 1/^^ integrated to zero. Putting these together gives the asymptotic behaviour 
of the boson mass: 

(53) 

which is exactly as predicted by the virtual ionisation arguments. 

We now evaluate the free energy. With the action in Gaussian form the functional integral can be performed exactly 
to give, 

= -\nJv[A*,A] e^" = - InDet"! (-/i^ + |^ - i^.) , (54) 

= Tr In (^-fib + - ii^Sj ■ (55) 

For given /i and U Gaussian fluctuations do not alter Tc, but this is not to say that they do not have any effect at all. 
Following Randeria et al. [ p9[ , who examined an s-wave continuum model, we calculate the fluctuation contribution 
to the fermion number density, n{fj,,Tc). For fixed n, the changes in fi{n,Tc) effectively renormalises Tc(ri). 

The fermion number density, n, is derived from the free energy, T = J-q +J^b, by differentiation with respect to the 
chemical potential, fi: 

n = riQ + On = (5o) 

The contribution from the free fermion action, J-'q, is 

no = l-^tanh(/3^/2)iV(e). (57) 

e 

In the BCS regime, |^| < 4, Sn is negligible and the mean field results are unaltered. As we approach the band edge, 
\fj,\ — 4, where our theory is not strictly valid, the fluctuation contribution starts to become important. In the Bose 
limit, l/^l > 4, the fluctuation contribution to n is considerable. In fact, in the U ^ oo limit n(/i,Tc) = 5n(/x,Tc). The 
boson contribution is, 

■^"^ ^^^^ -Aifc + q-q/2Af-io.,' ^^^^ 

= §^5]n6(-Mb + q-q/2A/), (59) 
q 

where the Bose occupation number is nb{x) = [e^^ — For simplicity we have ignored small contributions from 

terms proportional to the derivatives of d and m which tend to zero as U increases. 
The partial derivative of r is, 

which, for T=Tc, leads to the useful result that, 

lim ^ = 2. (61) 

Combining this result with Eq. ( |59| ) tells us that, in the low density, strong interaction limit, the fermion density at 
the critical temperature is, 

n(//,T,) = 2^nfc(q.q/2M), (62) 
q 

i.e., all the fermions are bound into Bose pairs of mass M . 

Unfortunately, there is a fly in the ointment: in two dimensions the integral in Eq. ( |6^ ) has an infrared logarithmic 
divergence due to soft q = modes . This problem can be circumvented be remembering that the cuprate materials 
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are not really two dimensional; rather they are highly anisotropic 3D materials. To take account of this we make the 
substitution 



2M 2M 2M 2M 2M 2M' ^ ' 

where = M/M^ ^ 1, reflecting the low mobility of electrons between Cu02 planes. With this substitution n is 

2 2 

always finite. Assuming Tc < we can safely take the limits of the q integral to infinity, giving, 



C(3/2) / 



(64) 



where (,{■) is the Riemann Zeta function. Aside from the 1/a prefactor which accounts for the mass anisotropy, this 
is exactly the result we would expect for a 3D gas of mass M bosons with Bose-Einstein condensation temperature 
Tc. Rearranging this equation gives an expression for Tc at fixed density: 

/ x2/3 

Tc = 2.09^^^. (65) 

From this equation we can see the importance of the mass: a large mass means a low Tc and vice versa. Thus, for 
the large U limit in which we are interested, increasing the coupling, which would "normally" (i.e., in mean field talk) 
increase the critical temperature, has the paradoxical effect of lowering Tc due to the localising effect it has on the 
pairs. 

We expect to see analogous behaviour in our extended Hubbard model as we increase the coupling, but instead of 
becoming localised on a site, our bosons should become localised on bonds, becoming more and more massive, i.e., 
less mobile, as V increases. 



B. Nearest neighbour (non-local) pairing 

The analysis for the extended Hubbard model proceeds along similar lines as the on-site model but is complicated 
by the existence of two species of bosons. As and Ad, or, equivalently, A^ and Ay. In the large V limit we find, 

(66) 



(67) 

For nearest neighbour hopping, these give the following bosonic masses: 

{Mlp,Mlp)^-^{ (1,1) a/3 = dd, (68) 
^ y (1,-1) a/3 = sd. 

The anisotropy between the s and d boson masses arises because the form factors for the s- and d-wave bosons, 
and C^C*: differ by the sign on cos(fc2;) cos(fcy). (This term does not enter the sd cross term.) This is best understood 
in terms of the x and y fields. In that basis we see that the term comes from the process, shown in Fig. ^(i), by which 
an X field becomes a y field. Obviously such a process is sensitive to the relative phases of the two bonds: zero phase 
difference is optimal for the hopping to occur so the s-wave phase is preferred. 

At this point, a pertinent question is, what happens if we change the bandstructure? As an example, consider 
adding a next nearest neighbour hopping, t', to the Hamiltonian so that the band energy is, 

Ek = —2t {cos{kx) + cos(fcy)) — At' cos{kx) cos{ky). (69) 

The diagonal hopping allows the x ^ y process to occur without a potential energy cost: there is no intermediate 
state where V is lost. See Fig. ||(ii) 



V^oo V 



V ^C"(k)C''(k) 1 



lim = \ " 

y^2TO^^"l6iV^ m*,(k) il' 

V ^C"(k)C''(k) (I , 2ek 



16iV^ <(k) 



,2 ' ,,3 
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In terms of the mathematics, the t' hopping means that the integrand for the 1/ fJ.^ term in the expansion of 1/^^ 
goes hke ztt' cos^ (k^) cos^ (ky) which has the same sign over the whole range of integration and is therefore finite 
(c.f. the nearest neighbour hopping contribution at this order which has only odd powers of cos(fc2:) and cos(fcj,) and 
therefore integrates to zero). The masses calculated with this term are: 

1 _ / -i(9t2 + t'V) a^s, 



If the product t'V is much smaller than the nearest neighbour hopping can be safely ignored. However, if this is 
not the case then the corrections can be of great significance. For instance, if t'V ~ t"^ then the d-wave solution is 
unstable (has negative mass), or if t'V ~ —9t^ then the s-wave solution is unstable. In the limit \t'V\ ^ t^ only one 
flavour of boson is stable and its mass is independent of V: 

lim M(i,<') = 4- (71) 

Including a sufficiently large diagonal hopping, t' , obviates the need for pair hopping to occur via virtual ionisation. 
Once this requirement has been lifted, the lattice model behaves exactly as we would expect a continuum model to 
behave (as shown in Fig. p. 

Assuming for the moment that both boson modes are stable (i.e., t' = or small) we resume our calculation. The 
Gaussian free energy is, 

= Tr In (^-^f + ^ - tiu}j , (72) 
= ^ln(A+A„), (73) 

where we have used the cyclic nature of the trace to diagonalise the action with respect to its flavour indices. The 
eigenvalues, A±, are the energies of the eigenmodes of the system which in general are linear combinations of s and d 
(or X and y) modes. In the large V limit, the eigenvalues are, 

A± = -i(A^L+MS.) + ^Jg^--., (74) 
where the masses, which depend on the angle 9 = tan~^(^), are, 



M± (9) M 



5±A/l7cos4(6')+3Ocos2(0)sin^(6') + 17sin'*(6l) , (75) 



with M = Mdd = V/t"^ as a convenient mass unit. 
With these definitions the bosonic number density is, 



n{^Ji,T,) = ^ p d9[n+{9)+n^{9)], (76) 

where the densities of the two modes, as functions of angle in q-space, are 

^±(^) = A/^^(M±(^)T.)t (77) 
V TT a 

We see from this that for given n and Tc, the mode with the largest mass will be the most populous. This is graphically 
illustrated in Fig. || which shows the relative contributions of the two masses to the angle-dependent density. 

The flavour of boson with the smaller of the two masses, M^{9), contributes negligibly (~ 3%) to the total number 
density, therefore the physics is determined by the boson with the largest mass, M_(0). From Fig. ^ we see that the 
M_ (9) bosons are not distributed evenly throughout space - they prefer to lie along one or other of the coordinate axes. 
Relative to an isotropic state, the bosons have been shifted away from the diagonals and onto the axes. For M_{9) 
as derived here, the effect is not very large (~ 10% of the total), but by adjusting slightly the masses, for instance by 
setting Mss — Mdd, nearly all the weight can be shifted to the axes, leading to a very pronounced anisotropy. 
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V. DISCUSSION 



To understand how this might relate to the cuprates it is helpful to take a step back. In the derivations above we 
have assumed that all the electrons are bound into Bose pairs, which is reasonable provided Tc/V <C 1. However, 
for higher temperatures, where T/y ~ 1, this assumption is not valid as thermal excitations lead to the break-up of 
pairs. In the extreme high temperature regime almost all the electrons are free. Now, starting in the high temperature 
regime, as the temperature is lowered, the pairing fluctuations become more important and bosons start to form. In 
the Gaussian approximation, the onset of fluctuations is gradual, with a maximum at Tc. T* can be interpreted 
as the temperature where the fluctuations become significant, which is not necessarily the same as T*^^. As the 
Bose number, 5n, increases, the fermionic spectral weight is correspondingly decreased according to Eq. ([5^). Given 
the strongly anisotropic Bose number in Fig. ^ we can expect a similar anisotropy in the pseudogap. Interestingly, 
the Bose number is maximum along the axes and minimum on the diagonals which is just the same as the angle 
dependence of the pseudogap observed by ARPES measurements. Lowering the temperature further leads to the 
creation of more bosons and the pseudogap becomes deeper. Finally, at the Bose-Einstein condensation temperature 
a real d-wave superconducting gap will start to open, the maxima of which "accidentally" coincide with the pseudogap 
regions where there are most bosons. 

It is important here to make two comments. Firstly, the pseudogap which we observe is not a superconducting 
effect in the sense that its origin does not lie in the spontaneous symmetry breaking of the superconducting order 
parameter, but rather in the formation of two-particle states which "eat" the available single-particle spectral weight 
in certain directions. Secondly, although the angular variation of the pseudogap is intimately related to the presence 
of d-wave pairs, the relation between the two is not as simple as might first appear. A single flavour order parameter, 
either s-wave or c?-wave, must have an isotropic mass tensor (for d-wave this is shown as the dotted circle in Fig. 
^ , otherwise the theory will not be invariant under the square symmetry of the lattice. The angular anisotropy of 
the pseudogap comes not from the symmetry of a single order parameter, but from a coupling of the d-wave to the 
s-wave bosons via the anisotropic gradient, A^(V^ — V^)As -I- c.c, which has the full square symmetry of the lattice. 
Equivalently we could say that the angular anisotropy is a reflection of the anisotropy in the masses of the a; and y 
bosons: an x boson prefers to travel in the x direction and a y boson prefers to travel in the y direction. Were the x 
boson to be totally confined to the cc-axis, and the y boson to the y-axis, in which case Mss = Mdd = M , the angular 
distribution of the bosons would be sharply peaked along the axes and zero on the diagonals. 

We have seen that it is necessary to have both s- and d-wave order parameters if the pseudogap is to have any 
angular dependence, so it is natural to enquire about the nature of the superconducting state formed below Tc. In 
the preceeding calculations it has been assumed that in the low density, large coupling limit, the critical temperatures 
for the s- and d-wave fields become degenerate. Now, just because two order parameters have the same Tc does not 
mean that they will be present in the condensate in equal measure. Ignoring the gradient coupling for the moment, by 



analogy with Eq. (77) it is easy to see that, in a mixture of the two pure states, the heavier d bosons will dominate by 



3/2 3/2 

a ratio of Ms : — 1 : 27. Including the gradient coupling terms the modes are admixtures of s and d bosons. 

The relative components can be found from the eigenvectors of the flavour matrix: 

where the s component is, 

4± x/l7cos4(e') +30cos2(6i) ain^ie) + 17sin''(0) 

«?(^) = ,,,, . 2,,, • (79) 

cos^(0) — sm (t^j 

For the dominant mode. Fig. |l^ shows the variation of the s component as a function of angle. (For the subdominant 
mode the same graph is produced by plotting l/e+(0).) For momenta lying along, or near, the diagonals, the dominant 
mode has an entirely d-wave character. As we move away from the diagonals and towards the x and y axes, the s-wave 
fraction grows, but even at its maximum value it is an order of magnitude smaller than the d-wave. We conclude that 
the dominant mode is essentially d-wave, whilst the subdominant one is essentially s-wave. 

Thus, by including both the s- and d-wave order parameters in the same theory, as well as the gradient couplings 
between them, we have been able to produce a pseudogap state with a non-trivial angular dependence, qualitatively 
similar to that of the d-wave superconducting state. Although the angle dependent pseudogap state requires an s — d 
admixture, the superconducting state that it eventually gives way to has, for most intents and purposes, pure d-wave 
symmetry. 
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VI. CONCLUSION 



In this paper we have considered an extended Hubbard model with nearest neighbour attractive interaction as 
a phcnomcnological model for the high Tc superconducting cuprates. By introducing a tensorial version of the 
Hubbard- Stratonovic transformation we were able to derive an effective Landau- Ginzburg- Wilson action describing 
both extended s- and rf-wave superconductivity. 

The saddle point (mean field) approximation to the full action gave the phase diagram of the model as a function of 
chemical potential (filling) and temperature. Near the centre of the band d-wave solutions were found to be dominant 
whilst the extended s-wave solutions were dominant near the band edges. The two flavours of superconductivity were 
separated from each other by two second order phase transitions, in between which there was found to exist a mixed 
time-reversal symmetry breaking s ± id phase. The alternative s ± d phase was ruled out except in the presence of 
orthorhombic distortions. At large values of the coupling, non-zero critical temperatures were observed beyond the 
band edges, signaling the creation of bosonic bound states. 

Using the Gaussian approximation to the LGW action, in the low density limit, where quartic terms representing 
interactions between bosons ought to be unimportant, we first reproduced the known large U behaviour of the on-site 
attractive Hubbard model, in which the bosons acquire a mass proportional to the coupling, U, due to the process 
of virtual ionisation. In the equivalent, large V, limit the extended model showed qualitatively similar behaviour 
provided only nearest neighbour hopping was considered. Introducing a next nearest neighbour hopping, i', allowed 
diagonal hopping processes, circumventing the need for virtual ionisation. It was shown that, in the limit \t'V\ 3> t^, 
there is only one flavour of stable boson, s-wave if > or d-wave if t' < 0, and that its mass is independent of V, 
just like in a continuum model. It was then argued that, for \t'V\ <IC t"^, the existence of two different bosonic modes, 
with anisotopic gradient couplings, could lead to an angle dependent distribution of bosons. This was then given 
a possible interpretation as the cause of the pseudogap state; two-particle states "eat" the single-particle spectral 
weight in certain areas of momentum space. Finally it was shown that, despite the mixing of the two different boson 
flavours, the superconducting ground state would be of predominantly d-wave character. 
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BCS U ^ Bose 

FIG. 1. Crossover from BCS to Bose pairing for continuum and lattice models. The full line shows the BCS mean field Tc 

evaluated for all strengths of the coupling, U. The broken lines show the effects of fluctuations on Tc. 



o 




o 










s- wave 


d- 


wave 



15 



(i) 



FIG. 2. s- and d-wave order parameters. 
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FIG. 3. (i) Quadratic and (ii) quartic terms in the expansion of Trln (1 — VGo). 
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FIG. 4. Weighted densities of states (WDOS) above Tc for the normal state and the extended s- and d-wave states. 
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FIG. 5. Mean field critical temperature, Tc, versus chemical potential, /i, for s- and i 
(b) V = 5t and (c) V = 8t. 
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FIG. 6. Plot of the stable phases of the Landau free energy as a function of the two dimensionless parameters, Us and 
The dashed line shows parameters calculated from extended Hubbard model parameters near the tetracritical point. 
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FIG. 7. (a) Real and (b) imaginary parts of tlie frequency coefficient, d{Tc,n), for extended s- and d-wave superconductivity 
atV = 16t. 
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